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Abstract 

The application of the locally harmonic approximation to large scale collective motion is 
briefly reviewed. Particular emphasis is paid to issues which might be useful in the more 
general context, or which are specific to our treatment, like there are: Self-consistency 
between collective and intrinsic motion, a T dependent coupling constant, the inclu- 
sion of quantum effects for the dynamics of fluctuations. Finally, open problems are 
addressed, like the use of thermodynamic concepts at smaller excitations, the ergodic- 
ity question, the variation of transport coefficients with the nuclear shape due to shell 
structure, the limits of the LHA at smaller T. 

1 Introduction 

In this talk I would like to briefly report on the transport theory which has been worked out on 
the basis of a suitable application of linear response theory. In the second part I shall address 
some of the most stringent and largely unresolved problems one generally faces in dealing 
with nuclear transport at smaller excitations. Details will be presented only to outline the 
general concepts. For more extensive studies I like to refer to my recent review article as 
well as to the original papers, especially to those of the more recent years which have been 
published together with F.A. Ivanyuk, D. Kiderlen and S. Yamaji and others. 

*Prepared for the RIKEN Symposium on "Dynamics in Hot Nuclei", Tokyo, March 1998 
^e-mail: hhofmann@. physik.tu-muenchen.de 

*http:// www . phy sik . t u- muenchen . de/tumphy/d/T36/hofmann.html 
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1.1 The basic approximation scheme 

Unfortunately, still too often theories of nuclear collective motion are considered separately 
from transport models. This is unjust both in the light of present day experiments as well 
as from the more general point of view. Indeed, it is probably only for rotational motion 
exactly along the Yrast line that intrinsic excitations of the nuclear system do not play any 
role. Moreover, one should not forget that one of the first theoretical descriptions of the 
prime example of large scale collective motion, nuclear fission, was given by Kramers as 
early as 1940 and which was based on the picture of "transport" in collective phase space. In 
his famous equation there appear such concepts as dissipation and stochasticity. (The reader 
who is not familiar with this equation is asked to wait for a short while; we will come back to 
it later). Kramers applied this equation to the decay rate for fission, to find a generalization 
of the famous Bohr- Wheeler formula. In these days the origin of dissipation was attributed 
to the strong " correlations" among the nucleons, as they can be understood within or follow 
from N. Bohr's compound nucleus. Below we shall try to elaborate on such a point of view 
on the basis of our present day understanding of nuclear dynamics. First we want to explain 
why and in which way large scale motion may be described by a linear response approach. 

The basic element is the observation that motion in the collective phase space may be 
expressed in terms of propagators by writing for the time evolution of the density distribution 

/(Q, P,t) = J dQodPo K{Q, P, t- Qo, Po, to) fiQo, Po, to) (1) 

with limt^t, K{Q, P, t; Qo, ^o, ^o) = 5(Q - QitoMP - P(to)) 




Q 

Figure 1: Time development of global motion within the propagator method. For the time 
step from t2 to ^3 the related Langevin method is indicated. 

Here K{Q, P, t; Qq, Pq, to) may be interpreted as the conditional probability for the system 
to move from Qo, Pq at to to Q, P at time t. On both sides of this relation the distribution /, 
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the "joint probability", may be replaced by conditional probabilities defining the transition 
say from a to to the final time t through an intermediate step at ti. The resulting relation 
is nothing else but the Chapman- Kolmogorov equation. (For a discussion of such general 
properties we may refer to the book by van Kampen 0). 




AQ On Q 

Figure 2: The LHA to a typical fission potential, which may either be given by the internal 
or the free energy. 

The procedure just described may be repeated introducing small time steps 6t. They only 
must be chosen longer than the typical microscopic time r below which an averaging over the 
microscopic dynamics of the nucleons would become meaningless. A linear response approach 
becomes possible whenever the 6t is small on the collective time scale. Consequently, the 
K{Q, P,t] Qq, PQjto) may be constructed within a locally harmonic approximation (LHA), to 
collective motion. If this propagator spreads over a limited regime of phase space only, within 
the time lap St, a Gaussian form will do. A pictorial view of this method is given in Figs.l 
and 2. 

1.2 The equation of motion for the propagator 

To approximate the propagator by a Gaussian means to write: 

KG{Q,P,t-Qo,Po,to) = (27r)-M-V2 exp{-^Y.i^^-^"ii))i^k-X^it))A^kit)} (2) 

i,k 

with the A(t) being the determinant A{t) = det(^^^(t)). The vector X defines a point in 
collective phase space, and the Kq represents a distribution, centered at the trajectory X'^(t), 
which starts at Pq,Qq, 

= (oil) = ({«)!) = 
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and which fulfills the following set of equations: 



Here we have introduced the variable q = Q — Qm- The Qm defines the position of the 
local oscillator, more precisely the position of the extremum of the local, effective collective 
potential. It is a linear function of Q{to), and has to be defined by the condition that the 
constant term in the force disappears see Fig. 2 and below. 
The matrix Aik{t) is given by the widths of the distribution 

with the Eqq(^) etc. being defined as Ej^) = HQ - Q^{t)f)t = {Q^)t - (<5'(^))^- These 
second moments fulfill the following set of equations: 

|E,^t) - ^Ejt) = (6) 
- j^EpJ^t) + CEq^t) + ^EqJ^t) = D,, (7) 

f^Epjit) + 2CE Jt) + 2^1:^^^) = (8) 

The inhomogeneities represent the diffusion coefficients. The initial condition J2q(lto) = 
J^qjito) = J^pjito) = warrants the Kq to start from the sharp distribution as required 

by®. 

This Kg satisfies the following transport equation 



d P d ^ d P r. d 

+ -TTFiCq + 7T7^— 7 + Dgp—— + D 



dqM dP ^ dPM ' '^"dqdP ^"^dPdP 



fiq,P,t) (9) 



This is easily verified by mere differentiation. The form of this equation resembles the one of 
Kramers with the exception that in his case the diffusion coefficients were given by Dpp = 7T 
and Dqp = 0. Moreover, as he understood his equation to describe large scale motion, the 
conservative force was given by the derivative of the full potential, —dV{q)/dq, rather than 
by -Cq. 



2 The derivation of the EOM within the LHA 

In this chapter we are going to describe briefly how we may obtain the transport coefficients 
as they appear in (H), or in the equations for the flrst @ and second moments (^H). To 
this end let us assume to be given a Q-dependent Hamiltonian of the type H{xi,Pi,Q) = 
HspiQ) + Vj'^J{xi,pi) The term Hsp{Q) represents independent particle motion in a deformed 
shell model, with Strutinsky renormalization meant to be included. The V}^^{xi^pi) stands 
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for a residual interactions, considered finally by way of a complex self-energies of the single 
particles: S{uj ± ie,T) = S'(co',T) =p |r(u;,T). Microscopic expressions are then calculated 
after replacing the single particle strength gk{uj) = 2n 6{huj — e^) by 

with the fi being the chemical potential. The F is understood to represent " coUisional damp- 
ing" . In numerical computations, the following values have mostly been used for the param- 
eters entering here: Fq = 33 MeV and c = 20 MeV. 

2.1 Average motion 

2.1.1 Secular equation from energy conservation 

An equation of average motion may be obtained from energy conservation f^, regarding 
the nucleus as an isolated system. One may thus write = ^-Etot = (5<^^^^^^^ 
Q (P{xi,pi,Q)^ ^, to get the equation of motion for Q{t) after expressing the average 
(^F{xi,pi, Q))^ as a functional of Q{t). Following the scheme of the LHA one may expand the 
H{Q) around any given Qq to have: 

H{Q{t)) = H{Qo) + m) - Qo)F + -m) - Q^f ( ^(Qo)) (11) 

The effects of the coupling term {Q{t) —Qq)F may then be treated by linear response theory, 
exploiting as a powerful tool the causal response function x 

^^t - s) = Q{t - s)'-tT {p^,{Qo,To)[F\t),F\s)])=2zQ{t-s)x"{t-s) (12) 

Here, the time evolution in F\t) as well as the density operator are determined by H{Qq). 
The is meant to represent thermal equilibrium at Qq with excitation being parameterized 
by temperature Tq or by entropy Sq. 

The equation of motion may finally be brought to the form k~^q{t) + x(t — s)g(s)ds = 
whose Fourier transform leads to the well known form of the secular equation for the possible 
local frequencies of the harmonic motion: k~^ + x{^) = 0. Its solution will be discussed below. 
The q{t) is given by 



q{t)=Q{t)-Q^ with ^^(^'^o) 



dQ 



^ d'E{Q,So) 



Qo 



(Qm - Qo) = (13) 



and Qm represents the position of the center of the local oscillator and the k a coupling 
constant, which for the case of zero temperature is known from the Copenhagen version of 
describing harmonic vibrations 0. At finite thermal excitations things become more com- 
plicated, where besides the shape variable itself also the parameter specifying the state of 
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the system may change in time. Choosing entropy, however, one may argue this change to 
happen only in second order in the velocity. For harmonic motion this is beyond the order to 
be considered. Keeping 5" constant one then obtains M 



+ X(0) (14) 

Qo 



The = " { ^dQ(^ ^ s adiabatic susceptibility and E{Q, Sq) is the internal energy. As 

the latter is somewhat difficult to evaluate one may introduce the free energy to rewrite (|1^) 
as 



+ x(0)-|^fJST) 1 (15) 



dTdQ 



Qo,To 



Notice, please, that for constant entropy S one may, for the quasi-static case, relate the 
temperature change to the one of the collective coordinate itself by applying general rules 
from thermodynamics: dT = {dT/dQ)gdQ. 

Another point which we want to elaborate further below concerns the difference between 
the adiabatic susceptibility and the static response, which can be shown to be either positive 
or zero. It vanishes at zero thermal excitations but at finite temperature it does so only for 
ergodic systems, a problem to which we shall return below: 



X^" - X(0) = (16) 



2.1.2 The collective response 



Let us introduce a coupling qextF" to a time dependent "external field" fextit) by adding 
fext{t)F to our Hamiltonian: H' = H{xi,pi,Q) + feyct(t)F but where H{Q) is taken in the 



approximate form (|TTD . The coupling is chosen to have the same form as the one between the 
two "subsystems" of collective and nucleonic degrees of freedom, with only the {Q(t) — Qq) 
replaced by the "external field" /ext- The "collective" response function may now be defined 
as: 

S{F)^ = -Xcoii(^)/cxt(^) (17) 

where for convenience we have switched to frequency representation. The difference between 
the S{F)^ and the {F)ij should at best have a term proportional to S^u), representing some 
static force in the time dependent picture. The phrase "collective" response function is chosen 
because i) this function now includes the collective excitations as well and ii) we will later on 
consider the latter as the only interesting ones. 

The construction of the Xcoii(i^) can be done following the derivation of the polarizability 
function for electric media as introduced by Clausius and Mosotti, and as described in ^ for 
the nuclear case at zero excitation. We need three equations: Besides the basic definition ([T7|), 
we have to have an equation which relates {F)t to {Q(t) — Qq) and we need to express {F)t 
by the total perturbation consisting of both external and induced fields. We want to perform 
this construction assuming ergodicity for our nucleonic system, i.e. requiring the condition 
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( ITBD to be fulfilled. Exploiting ([I^ the relation between the average {F)t and the deviation 
of Q from Qq becomes 

k{F)t = Q-Qo (18) 
Finally, the collective response function defined by (|l^) takes on the form 

Xcoii(c.) = -r^Pf-, (19) 

The relation (p!8|) reminds one of a self-consistency condition, like one would find treating 
an effective two-body interaction {k/2)FF in mean field approximation. However, such a 
notion is justified on more general grounds. It is this relation (|l^) which expresses most 
clearly that the correct interpretation of the Q{t) must be to consider this quantity an internal 
variable — measuring some properties of the nucleus and clearly to be distinguished from truly 
external fields like the /ext- 



2.1.3 Transport Coefficients and the distribution of collective strength 

The procedure of defining the transport coefficient of inertia, friction and local stiffness is 
as follows: For selected peaks of the strength distribution one fits the dissipative part of the 
oscillator response Xosc(^) = (Mcu^ + i'juj — C) ^ to the Xcoii(i^) calculated microscopically. 
If one is interested in the transport coefficients for the g— motion one first transforms to the 
response function Xqqil^) i^i fhe g-mode, defined through Xqqi!^) = '5g(a;)/5(— gext(i^) |ijoxt=o- 
This can easily be done with the help of the linear relation: kS{F)t = SQ{t) = q{t) which 
implies Xqqi^) = ^^Xcoii(^)- Summarizing, this fit then implies the substitution (with cuf = 
±£ - iT/2) 

Xqq{uj) =^ Xosc(^) = JTT-f — + (20) 

M[UJi — uJi ) \UJ — uJi uj — uJi J 

in a certain range frequencies which encompasses the peak in the strength distribution whose 
associated mode ui one wants to treat explicitly. Notice, please, that in this way we actually 
solve the secular equation for arbitrary damping. In the case of fission it is believed that it 
is mainly the low frequency behavior which matters, for which reason it is the lowest peak 
which is treated in this way. For an application to giant resonances see |^ for the case of the 
quadrupole and the talk by S. Yamaji at this meeting for monopole vibrations. 

Within the formulation presented above the strength distribution may change dramatically 
with temperature 0, p[, 0. This is mainly due to the temperature dependence of the 
coupling constant, as determined by the stiffness C(0) = d'^E{Q, So)/dQ^ of the static energy 
(see (0)). For the quadrupole as well as for the local motion along the fission path, the 
C(0) may well decrease by more than an order of magnitude when T is increased from 
to about 2 ■ • • 3 MeV. As a consequence, at the higher temperatures practically all strength 
concentrates in one dominant low frequency mode, with the associated transport coefficients 
reflecting strongly over- damped motion. This feature has been demonstrated in (see Fig..l 
of this reference or Fig. 5. 1.1 of [IH). As may be expected from such a strength distribution. 
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the inertia of the low frequency mode decreases from the typical cranking model value at low 
excitations to the one given by the energy weighted sum rule at large T 0. 
Conveniently the transport coefficients are expressed by the following ratios 

'CI 7 



M 



M 



V 



2./M I C 



(21) 



with which the frequencies the two frequencies may be expressed as cu^ = ±£ — 



w (^±-v/signC — 7]'^ — ir]j. Numerical values from microscopic computations can be found in 
[§]) m ^'^'^ iini- A typical calculation of F is shown in Fig. 3 (from ||10|)- If evaluated at the 
potential minimum and at the barrier, typically the hw is of the order of 1 MeV, indepen- 
dently of T; conversely the rj increases from about 1 at T = 1 MeV to about 3 at 5 MeV 

0- 
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Figure 3: The inverse relaxation time F 
temperature (indicated in the Figure). 



7(a;i)/M(a;i) as function of deformation and 



The behavior found in our model for the strength distribution may be understood in the 
sense of reaching a macroscopic limit with increasing T. Within such a limit there is no room 
any more neither for quantal coherence in the sense of collective motion nor for the typical 
single particle excitations. The nucleus just behaves like a drop of a strongly viscous nuclear 
fluid (see also [|l^ , [|12] ) . This behavior is in clear contrast to the one found for common RPA 
at finite T, for instance like in the so called "thermo field theory" (see talks at this meeting 



andc.f. [|T3l). 



2.1.4 The equation of average motion and the transfer of heat 

After the substitution (^) the equation for the g-mode attains the form: 

M(cui) q{t) + 7(cui) q{t) + C{uJi) q{t) = -gext(t) 



(22) 
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It is of differential nature although we have not really made use of any Markovian approx- 
imation. All we did is to restrict ourselves to just one mode assuming that the associated 
peak can be approximated by a Lorentzian. For these reasons we are actually dealing with 
Non-Markovian transport coefficients representing dynamics in a whole range of frequencies. 
They are determined uniquely for under- as well as for over-damped motion, and for positive 
or negative stiffness. 

Let us turn to the energy balance once more. One finds the equation 

d d (Me Cq'\ .2 d dS 

no matter whether the external field is "switched on or off' For (damped) self- sustained 
motion the equation dEtot/dt = may thus be rewritten as 

d _ d ( M{uJi) .2 C{uJi) A ..2_rj.d , 

which correctly expresses the exchange between collective motion into heat. 

2.2 The dynamics of fluctuations of the collective degrees of free- 
dom 

So far the collective variable Q was just a c-number classifying the deformation of the mean 
field. To be able studying collective fluctuations we first must raise Q to become a genuine 
dynamical variable. As we are interested in quantum effects this generalization has to be 
performed in operator sense. Simultaneously, we need to introduce the associated momentum 
n, and we must find a decent Hamiltonian for the total system. Such a Hamiltonian may then 
be used to derive the transport equation for the collective density operator rf, the Wigner 
function of which may be associated with the distribution in collective phase space mentioned 
before. Evidently, such a procedure ought to account for the basic and important features 
we have been able to consider for average motion, the most stringent one perhaps being self- 
consistency. If at all possible, this ambitious program can be accomplished only in a restrictive 
sense; here it will be the LHA approximation, once more, which is exploited. 

2.2.1 A suitable Hamiltonian for the total system 

We would like to have a Hamiltonian which splits into parts of nucleonic and collective motion 
plus a coupling between both subsystems: H = ifnuci+-^coup+-^coii- This goal may be achieved 
by applying the Bohm-Pines method, which originally has been invented to treat collective 
modes in the electron gas. One obtains: 

n = ^(x„p„ Qo) +kIlF-^iQ- Qo)F +I^ + ^(2f3+k)iQ- Qof (25) 
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The unperturbed collective inertia is given by the sum rule value, namely 

— =te{ [P, F] )q„ P = t [HiQo), F] (26) 

mo 

The n is the canonical momentum satisfying the commutation rule [Q, fl] = ih. 

The intrinsic part obtains a most convenient form: It is the same H{xi,Pi,QQ) which 
appeared already in (|TT]) as the unperturbed part. The important point is that there is no 
reduction of the number of nucleonic degrees of freedom. This implies that we will be able to 
work with the same nucleonic response functions as before. Of course, in order not to have 
too many degrees of freedom there will have to be a subsidiary condition, which turns out to 
be 

kF-{Q- Qo) = (27) 

For average motion it reduces to the form (|I8|) found before. Notice the two coupling terms: 
The first one is of similar structure as that in ([TT|). In addition there is term involving P, 

which is multiplied by the same factor F which appears also in the unperturbed inertia. It 
is this intimate relation between the coupling terms and the unperturbed collective Hamilto- 
nian which in the end assures self-consistency, in the sense of the collective coordinate being 
related to the intrinsic field F. The alert reader will have noticed that there is the additional 
parameter /?. In the end the latter is chosen as {P + k)/k'^ = C = —MuiUi, Very naturally 
this goes along with the re-introduction of the kinetic momentum, such that for the averaged 
dynamics one regains the form d). 

For undamped motion the situation is like in common RPA: The Hamiltonian can be 
transformed in such a way that the coupling between collective and nucleonic degrees of 
freedom effectively shows up in renormalized transport coefficients for inertia and stiffness, 
after the intrinsic degrees of freedom have been averaged out. The corresponding frequency 
satisfies the secular equation of RPA which corresponds to the separable interaction {k/2)FF. 
This equation has the form given above but where in the x only the reactive part plays a 
role. For damped motion the situation is more complicated. To get information about the 
dynamics of collective fiuctuations one needs to invoke projection techniques, like the one of 
the Nakajima-Zwanzig method. 

2.2.2 A non-perturbative Nakajima-Zwanzig approach 

One starts with the von Neumann equation for the Hamiltonian (pSj) and defines the reduced 
density d by averaging out the intrinsic degrees of freedom. Its equation has a term which is 
non-local in time. The integral kernel /C contains an operator for the time evolution of the total 
system, which commonly is treated in second order perturbation theory, the argument being 
that JC does already contain the couphng twice. However, in order to meet the standard set 
up at the beginning of this subsection one needs to do better. The new scheme to be applied 
may be read of from the structure of the collective response function ([l9|) : On the one hand 
one may exploit low order perturbation theory with respect to the effects of the coupling on 
the intrinsic degrees of freedom, for which reason the nucleonic response function x appears. 
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Conversely, with respect to collective motion an infinite order is considered, as seen by the 
fact that the coupling constant appears in the denominator. To translate this scheme to the 
NZ approach essentially one replaces the unperturbed time evolution operator by an effective 
one: 

exp[-i(Lcoii + I/nucOs] =^ exp[-i(£[;^i{ + Lnucl)s] (28) 

Here, time evolution is formulated in terms of Liouvillians. Fortunately, the effective one, 
CfJJi, need not be constructed explicitly if we are interested only in the Gaussian form of the 
local propagators introduced in sect.l. It suffices to exploit the quantal fluctuation dissipation 
theorem, which can be proven to be satisfied within the modified Nakajima-Zwanzig approach. 

2.2.3 Diffusion coefficients and the equilibrium of the damped oscillator 

Let us turn our attention to equilibrium fluctuations first. Imagine the total system to be 
close to equilibrium at some given temperature T. The collective fluctuations may then be 
calculated by applying the quantal fluctuation dissipation theorem, plus a few general relations 
of response functions. For the three possible combinations of q and P one gets: 

^ 2 1 

^PP M^PP 99 ]\^^ip IP ^^gq M^PP ^ ^ 

It is quite easy to convince one-selves that these relations are in accord with (|]-§) for stationary 
situations. The equilibrium fluctuations are determined by 



= / ^ ^eoth xL(-) = M'I ^ ^-oth (30) 

For the kinetic momentum we have used the fact that P = Mq. The cross fluctuation vanishes 
because q and P behave differently under time reversal. We are now going to demonstrate first 
that for two limiting cases these expressions lead to results anticipated on general grounds. 

(i) Limit of high temperature: Let us begin addressing a situation for which we expect 
the Einstein relation proper to apply. Suppose the temperature to be sufficiently high such 
that in the integrals appearing in ( pOf ) we may effectively replace ^coth(^u;/2T) by 2T/uj. 
Then for the average kinetic and potential energy we just obtain the values known from the 
equipartition theorem, namely 

2M 2M 2 2^^ 2^«^ ^ ' 

For the diffusion coefficients given in (^) this indeed implies the forms given below eq.(^. 

(ii) Limit of zero damping: Simple results can also be obtained for the limit of zero 
friction. Take the oscillator response and let the imaginary part F of ujf go to zero. 
Then the dissipative part of the response function can be written as lim^^o Xosc(^) = 
{'k/2Mw) {5{ijJ — w) — 5{uj + w)). For such a mode, the average values of the kinetic and 
potential energies take on values given by the quantal version of the equipartition theorem, 
namely: 

- ^ ^-'i - ^ cothf— ) --(a^) - -V*^" r32) 
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Also, the diffusion coefficients then look like the ones obtained before with only T replaced 
by an effective temperature T*{w), namely 



qp 







Dpp = '^T*{w) 



with 



coth( ) 

2 ^2T' 



(33) 



This version of the generalized Einstein relation has been used in ref. flj] to describe the "fast" 
mode of charge equilibrization in heavy ion collisions. 

(in) General case: The integrals in (|30| ) may be evaluated in terms of infinite series 
by applying the residue theorem, once the hyperbolic cotangent has been developed into the 
uniformly convergent pole expansion involving the Matsubara frequencies. For evaluating Y^^p 
the integral must be regularized, for instance by introducing a frequency dependent friction 
coefficient which drops to zero beyond the so called Drude frequency ujd, for details please see 
Typical numerical results are shown in Fig. 4, where the "equilibrium fluctuation" YTqq 
is shown together with both diffusion coefficients, all quantities normalized in self-explaining 
fashion. They are plotted as function of the local stiffness C, which becomes negative when 
the local modes become unstable. It can be proven that this happens exactly when the 
(unperturbed) stiffness C(0) of the static energy turns negative. The figure demonstrates 
that our way of calculating the diffusion coefficients is continuously possible independent of 
the sign of the stiffness. It is only for smaller temperatures that one runs into difficulties for 
unstable modes; we will address this point later on. Incidentally, one also observes that for 
C < the YTqq itself is negative. This does no harm as for unstable modes this quantity itself 
has no physical meaning. 



o 




Figure 4: Diffusion coefficients (diagonal: thick fully drawn line, off diagonal: dashed line) 
and analytic continuation of the equilibrium fluctuations in q (thin full line) as functions of 
the stiffness C, for fixed values of M and 7. 
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3 Open questions 



In this section we want to address a few of the most stringent but unresolved problems; to 
some extent they have been mentioned already in papers cited above. 



3.1 Nuclear thermodynamics 

It has turned out unavoidable to benefit from the concepts of thermodynamics when classifying 
excitations of nuclear systems, knowing the latter not only to be small and isolated but to 
be unstable as well at already such small values of overall excitations as 10 MeV. Often, 
in the spirit of a shame-faced concealment, one uses (or perhaps misuses) the excuse that 
temperature (or the total energy) only parameterizes the level density. To me it is unclear 
what such a phrase really means, considering the fact that even in thermodynamics proper 
concepts like entropy and temperature only derive from certain interpretations of the level 
density, even if applied to truly macroscopic systems. Rather, one may apply the concepts of 
thermodynamics in literal sense, if one only adds the question about possible errors, in the 
quantitative sense. Certainly, such a point of view pre-supposes that the system cquilibrizes at 
all. This question simply depends on time scale: There must be sufficient time for most of the 
degrees of freedom available in a certain process to envisage " all" its possible configurations 
before the nucleus disintegrates in one way or other. Knowing that this condition may not 
always be easy to fulfill, in particular at the somewhat higher excitations, we want to start from 
such a hypothesis to concentrate on the problem about the use of the concept of temperature 
or the canonical distribution. Looking back at the various forms of the diffusion coefficients 
this question certainly appears to be relevant also in the more physical sense. 



3.1.1 Fluctuations of basic thermal quantities 

(i) The canonical distribution: Here one assumes the temperature to have a definite value, 
fixed by the internal energy £. The latter, however, is not determined precisely but exhibits 
a fluctuation which by way of the fundamental formula AS = T'^dS/dT = T'^C can be 
expressed through the specific heat C = TdS/dT = d£/dT. The entropy S is subject to the 
uncertainty AS = C. To get the right order of magnitude of these quantities it is totally 
legitimate to use the "Fermi gas" model for the level density: f2{S*) oc exp with 
S — 2y/a£* and £* — £ — Eq — aT'^ being the excitation energy. The specific heat becomes 
C = 2aT — S such that the relative fluctuations can be expressed as 

1 A£* AS 1 2 /2MeV\^/^ 

As indicated on the very right, a level density parameter has been used which is about twice 
as large as that for the genuine Fermi gas (free particles in a box). In this sense our estimate 
may be considered realistic for the nuclear case. From the typical values given below in the 
Table. 1 it is seen that the energy fluctuations take on considerable values! 
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(a) The micro- canonical distribution, of course, is exempt from this deficiency. It may be 
defined as p(e) = 1/G for E — A < e < E, with G measuring the number of levels inside the 
the interval considered and with p(e) = outside. As can be seen from Table. 1, even for a 
very small width A of the distribution there are so many states that a statistical treatment is 
unavoidable. Temperature may again be introduced through the relation = dS{e)/de Is, 
but now this quantity exhibits a finite fluctuation whose relative value is given by 



AT 



c 



1 A^* 
2 



(35) 



At smaller excitations this ratio may become uncomfortably large, a fact to which attention 
has been drawn already in [|15||. On the other hand, for an application to transport theory 
such a fluctuation does not hurt too much as long as the transport coefficients will not depend 
on temperature too sensitively. Such a situation may be given at T not smaller than 1 MeV 
Table. 1 

a) Dependence on A for ^* = 20 MeV, A = 0.01 MeV 



A 


50 


100 


150 


200 


T 


1.8 


1.3 


1.0 


0.9 


AS*/S* 


0.42 


0.36 


0.32 


0.30 


A/S* 


5 ■ 10-^ 


5 ■ 10"^ 


5 ■ 10-4 


5 ■ 10-4 


G{S*) 


3.7-105 


3.9 ■ 10^ 


4.8 ■ 10^2 


1.9 • 10^5 



b) Dependence on S* for A = 100, A = 0.001 MeV 



£* 


25 


75 


125 


175 


T 


1.4 


2.5 


3.2 


3.7 


AS*/S* 


0.34 


0.26 


0.22 


0.21 


A/S* 


4-10-5 


1.3 ■ 10-5 


8 ■ 10-6 


5.7-10-6 


G{S*) 


1.3 ■ 10^° 


7.5 ■ 10^0 


2.5 ■ 1028 


3.4 ■ 10^4 



3.1.2 The variation of temperature with deformation 

Applying transport equations to nuclear collective motion it is evident that temperature 
cannot be treated as being constant (see e.g. |T^, 0). One reason why T may change 
is because of the production of heat due to dissipation (as given by eq.(P^)). In practical 
applications this effect has commonly been taken into account (for an early treatment see 
[Q'lUll)- However, for an isolated system temperature may change even without dissipation. 
Within our formulation this effect is hidden in the way the coupling constant is calculated, 
and it is this feature which we are going to address now. 

In the hypothetical case of a constant temperature the coupling constant would be given 
by a form similar to (0), but with S replaced by T, the internal energy E replaced by the free 
energy and the adiabatic susceptibility by the isothermal one = ~ (d{F)'^^ /dQ^^ 
namely: 



T 



d'H.^ A , , T^_ d'J^iQ,To) 




(Qo) + (X(0) - X 



+ X(0) (36) 
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For such a process the quasi-static entropy would change as function of time. Therefore, if 
we beheve in our argument given before, this version (|36|) cannot be expected to correctly 
represent the physical situation. However, in many cases the difference is not very large. It 
can be expressed by the difference between the isothermal and adiabatic susceptibility which 
in turn may be calculated from derivatives of the free energy: 

Indeed, usually the term on the right hand side is not very big. Generally speaking, for 
T >~ 1.5 MeV one finds ^ - X^*^ < d^^{Q,To)/dQ\^j,^, perhaps except in the 
neighborhood of turning points of the potential energy. At small T, however, the assumption 
of T = const, may become questionable, safe for very special cases like vibrations about the 
ground state minimum of nuclei with doubly closed shells, whose deformation does not change 
with excitation. 



3.2 Ergodicity problem and the heat pole 

Above it had been mentioned that for ergodic systems the adiabatic susceptibility becomes 

,ad 



identical to the static response: x'^'^ = x(0)- Indeed, as shown in [^, (see also [|l|) it suf- 
fices to have (i) a non-degenerate spectrum of the intrinsic eigenstates as well as (ii) a narrow 
distribution in the occupation of these states. Evidently, both conditions will be difficult to 
fulfill within the pure single particle picture and the canonical ensemble. This statement may 
be made quantitative by invoking the correlation function il)"{u!) which relates to the dissipa- 
tive part x"{^) by way of the fluctuation dissipation theorem: hx"{uj) = tanh {huj/2T) 
As the derives from an expression similar to the one given in (|1^), but with the com- 

mutator replaced by an anti commutator, a 6 function type singularity appears at o; = 0: 



r27:6{uj) + Rr{uj)=T{x^ -x{0))27f6{uj) + R^i^) 



(38) 



with the Rip"{uj) being regular at u = 0. As written here, the pre-factor of this 6 function 
can be expressed by the difference of the isothermal susceptibility and the static response. In 
the sequel we like to call this contribution the "heat pole", with the ip^ being its residue. 
If evaluated within the pure independent particle picture one gets 



ipm 



J2 \Fkk\'^n{ek){l - n{ek)) = T J2 

k k 



dn{e) 
de 



=efe 



'djck - /i)' 



(39) 



In p] this 'ipfpra has been calculated as function of T. At low T it increases rapidly to reach a 
constant value above about T = 2 MeV. This value overshoots considerably the one it would 
have for an ergodic system, namely x^ ~ X^'^- Unfortunately, this situation does not change 
when the coUisional damping is taken into account. In this case the heat pole attains a finite 
width 

but the strength of the Lorentzian essentially remains unchanged, see P] and Jll. (It may 
be noted in passing that the width increases almost linearly in T, following the simple 
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rule ~ 2r(/i, T) ^ 2T, if considered over the large range of T = — 6 MeV). Two 
reasons may be responsible for this feature: (i) The calculation still is performed within the 
(grand) canonical ensemble; (ii) the choice of the self-energies in ( p!OD does actually not lift 
degeneracies. 

It seems to me that the questions raised here suggest an interesting and important problem 
for future studies, not just for our formulation of transport theory but in the more general 
context. The problem at stake touches two critical issues of nuclear physics: the use of both 
the independent particle model as well as of the canonical ensemble. We know that both 
models are not really applicable, and the evaluation of static susceptibilities allows one to 
put the question on a quantitative level. Indeed, might one not expect that a computation of 
these on the basis of the compound model and for the micro-canonical ensembles would show 
ergodicity, in the sense of (|T^)? 

To demonstrate the importance of such questions for transport properties let us look at 
nuclear friction. Often the latter may be approximated by the so called zero-frequency limit, 
for which one has: 

, , dx"(uj) ij"(uj = 0) 2^/ ^ , ,N 

7(0) = = ^^^^^ = rv + ^^'^ 



duo 



In the expression on the very right, the first term represents the contribution from the heat 
pole, the second one from the rest. Estimating — x(0) through ( |5SD this component of 
friction turns into the one found first by Ayick and Norenberg within the model of DDD 



for larger temperatures. In Fig. 5 it is shown by the full and dashed lines (see p]) as function 
of T. They correspond to values of the parameter c of (p!0D put equal to c = 20 MeV and 
c oo, respectively. The curve with the heavy squares corresponds to the contribution _r7(0) 
of the remaining part of the correlation function. The horizontal line represents wall friction. 
As compared to experimental evidence from nuclear fission, the friction force given by the (^) 
(with all terms included) appears to be too big. Furthermore, one expects nuclear friction 
to increase with excitation, rather than to decrease, also in the regime of somewhat higher 
T. From our discussion above the culprit for such a misbehavior is easily traced back to a 
violation of "ergodicity", in computations which base on the model of independent particles 
in a (grand) canonical ensemble. Therefore, it has been argued in and |l| to "re-install" the 
latter by neglecting in (^Ij) the contributions from ~ x(0)- For the "experimental facts" 
mentioned see e.g. and as well as the talk by G. Rudolf at this meeting in which a 



new interpretation of such data has been presented. It is only fair to add that at present no 
definite answer to such questions is possible yet. However, there can be little doubt that they 
may be considered as very important ones, not just for nuclear dissipation, but to understand 
nuclear collective motion in the more general context. 



3.3 Fluctuation of transport coefficients with shape at small exci- 
tations 

As seen from Fig. 3, strong variations with Q occur at the lower temperatures of 1 and 2 MeV. 
In principle, they show up because of shell effects. Unfortunately, in practical applications 
they appear to be grossly overestimated, simply because of unphysical irregularities in the 



independent particle model. It seems that even pairing correlations do not suffice |25|] to 



smooth out these oscillations. This problem may hint to the need of using truly adiabatic 
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Figure 5: Heat pole friction as function of temperature (see text) 



states. The latter would be produced after considering in an improved way effects of the 
residual interaction Vj-^J{xi,Pi). Apparently, this feature is very much related to the one 
encountered before when discussing susceptibilities. The Ve^l^j^P*) be expected not 
only to lift degeneracies but to repel states, not just at the spherical configuration but for 
deformed states as well. 

The puzzle raised here may be addressed from a different point of view: It is very doubtful 
that such small scale variations may at all be "seen" by the phase space distribution which 
moves across the collective landscape. This f{Q,P,t) itself will have some natural width, 
both in Q as well as in P. Therefore, one may well question any use transport coefficients 
which exhibit such extreme variations in the collective degrees of freedom. Obviously this 
problem is related to the one raised by the late V.M. Strutinsky (see e.g. [0). He questioned 
any introduction of collective degrees of freedom which is not accompanied by a suitable 
averaging procedure. Please notice that the variations we are talking about here are of a 
scale much smaller than those one would see in the typical gross shell behavior. Indeed, the 
ones shown in Fig. 3 for 7/M are related to those which in the potential energy would occur 
in between the few maxima and minima of a typical fission landscape (cf. e.g. [l^])- As for 
the transport coefficient of nuclear dissipation, the reader may be reminded of the physics 
of the wall formula. The latter represents a certain macroscopic limit, in which shell effects 
are washed out totally (cf. and This is the reason why no dependence on Q survives 
which might be related to shell effects of any kind; this point is elaborated on in [|TT] . 
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3.4 Breakdown of diffusion across a barrier at small excitations 

It has been mentioned already at the end of sect. 2 that diffusion coefficients may be calculated 
even for unstable modes, which are present at C(0) < 0. This is by far not trivial as for such 
modes the "equilibrium fluctuations" which appear in (|29|)-(|30|) loose their physical meaning. 
Such a calculation requires a delicate extension of linear response theory and the FDT to 
instabilities; for an comprehensive discussion see Unfortunately the procedure breaks 

down at some critical temperature Tc below which the diffusion coefficient Dpp would become 
negative. This may be inferred from Fig.6, together with (P^). The formal reason for this 




0.2 0.4 0.6 0.8 1.0 

T/G3 



Figure 6: Analytic continuations of the (diagonal) equilibrium fluctuations to the case of 
negative stiffness (dashed line for the coordinate, fully drawn lines for the momentum), as 
functions of temperature for different values of friction. The insert shows the momentum 
fluctuations for a larger range of the temperature and the same values of rj. 

problem may be analyzed in the limit of weak damping: For C ^ —C the coth in the T* 
introduced in (|33| ) simply turns into a cot! In this limit the is larger than the so called 
"cross over" temperature Tq found for "dissipative tunneling" by exactly a factor of 2. With 
increasing rj both Tq as well as T^. become smaller (with their ratio changing as well). For 



more details see 13 



As may be deduced from Fig.6 the problem occurs at temperatures of the order of T ^ 
0.5 MeV or less. Thus it is in this regime where the concept of the LHA may be said to 
fail, in regions where instabilities occur. However, as demonstrated in above this our 



method is capable of accounting for quantum corrections to Kramers' decay rate of fission, 
which without any doubt is a highly non-linear process. 
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3.5 The problems with the entrance phase of a heavy ion colUsion 



The theory discussed above was based on the assumption that it is possible to devide the 
system into "fast" and "slow" degrees of freedom. The former are supposed to exhibit genuine 
relaxation to a " local" equilibrium on top of which motion of the slow degrees of freedom may 
be described exploiting a quasi static description. Such favourable conditions are hardly 
fulfilled at the very early stages of heavy ion collisions, even if the energy per particle does 
not exceed some few Mev. The composite system starts from two separated fragments at 
zero intrinsic excitation. It will take a finite time before a common potential builts up with 
respect to which it may become meaningful again to speak of thermal excitations. During 
this phase the collective motion which is manifested in the distance between the two center of 
masses of the two (more or less still existing) fragments cannot be associated to low frequeny 



modes. This feature has already been discussed in [2S]. There it has been demonstrated 



that a large part of the excitation may be attributed to high frequency transitions in the 
nucleonic response, which nevertheless may lead to a kind of "pseudo-friction" for relative 
motion. In non-Markovian equations of motion have been suggested to dwell with this 
problem in the so called "diabatic dissipative dynamics" (DDD). It is conceivable that such 
effects could eventually be treated by generalizing the LHA suggested above in such a way 
that one concentrates on high frequency collective modes, rather than on the lowest one. 
After all, in some schematic way, the local response calculated within the DDD model looks 
very similar to strength distributions where all more detailled excitations are neglected other 
than just one prominant giant resonance, see [Q, A combination of this picture of 

average motion with our way of handling the dynamics of fluctuations might allow one to 
improve on the calculation of the diffusion coefficients. This seems to be highly needed as for 
the entrance phase the latter cannot be calculated within the high temperature limit. 

Acknowledgements: The author greatfuUy acknowledges support from the Deutsche 
Forschungsgemeinschaft. Furthermore, he likes to express his deep gratitude to his colleagues 
F.A. Ivanyuk, D. Kiderlen, A.G. Magner and S. Yamaji for their close collaboration over the 
past decade. 



References 

[1] H. Hofmann, Phys. Rep. 284 (4&5) (1997) 137-380 
[2] H.A. Kramers, Physica 7 (1940) 284 

[3] N.G. van Kampen: "Stochastic processes in physics and chemistry ", North-Holland, 
1981, Amsterdam 

[4] D. Kiderlen, H. Hofmann and F.A. Ivanyuk, Nucl.Phys. A550 (1992) 473 
[5] H. Hofmann, Phys. Lett. 61B (1976) 423 

[6] A. Bohr and B.R. Mottelson, "Nuclear Structure", Vol. II (Benjamin, London, 1975) 



19 



H. Hofmann, S. Yamaji and A.S. Jensen, Phys.Lett. B 286 (1992) 1 

H. Hofmann, F.A. Ivanyuk and S. Yamaji, Nucl. Phys. A 598 (1996) 187 

S. Yamaji, F.A. Ivanyuk and H. Hofmann, Nucl. Phys. A612 (1997) 1 

F.A. Ivanyuk, H. Hofmann, V.V. Pashkevich and S. Yamaji, Phys. Rev. C 55 (1997) 
1730; |nucl-th/9701032| 



H. Hofmann, F.A. Ivanyuk and A.G. Magner, |nucl-th/ 9709043 



A.G. Magner and H. Hofmann, to be published 

H. Hofmann and D. Kiderlen, to appear in J. Modern Physics E; |nucl-th/9703056| 



H. Hofmann, C. Gregoire, R. Lucas and C. Ngo, Z. f. Physik A293 (1979) 229 

H. Feshbach, Physics Today (Nov. 1987) 9 and IEEE J. Quant. Electronics, 24 (1988) 
1320 

H. Hofmann and P.J. Siemens, Nucl. Phys. A257 (1976) 165 
H. Hofmann and C. Ngo, Phys. Lett B 65 (1976) 97 

C. Ngo and H. Hofmann, Z. f. Physik A282 (1977) 83 

D. Kiderlen, Diploma Thesis, TU Munich (1990) 

R. Kubo, M. Toda and N. Hashitsume, "Statistical Physics II, Non-equilibrium Statistical 
Mechanics", Springer, 1985, N.Y. 

W. Brenig, "Statistical Theory of Heat, Nonequlibrium Phenomena", Springer, 1989, 
Berlin 

S. Ayik and W. Norenberg, Z. Phys. A309 (1982) 121; see also: W. Norenberg, in "Heavy 
Ion Reaction Theory", ed. by W.Q Shen, J.Y. Lin and L.X. Ge (World Scientific, Singa- 
pore, 1989) p. 1 

D. Hilscher, I.I. Gontchar and H. Rossner, Physics of Atomic Nuclei 57 (1994) 1187-1199 
D.J. Hofman, B.B. Back and P. Paul, Phys. Rev. C51 (1995) 2597 
F.A. Ivanyuk and H. Hofmann, to be published 
V. M. Strutinsky Nucl. Phys. A502 (1989) 67c 

H. Hofmann, G.-L. Ingold and M. Thoma, Phys. Lett.B 317 (1993) 489 

P.J. Johansen, P.J. Siemens, A.S. Jensen and H. Hofmann, Nucl. Phys. A288 (1977) 152 



20 



